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ABSTRACT 

We demonstrate that magnetically-collimated protostellar winds will sweep ambient material into thin, 
radiative, momentum-conserving shells whose features reproduce those commonly observed in bipolar 
molecular outflows. We find the typical position- velocity and mass- velocity relations to occur in outflows 
in a wide variety of ambient density distributions, regardless of the time histories of their driving winds. 

Subject headings: hydrodynamics — ISM: jets and outflows — stars: formation — stars: mass loss 



L INTRODUCTION 

As they form, young stars emit powerful winds whose 
mechanical luminosity amounts to a fair fraction of the 
stars' binding energy. These winds are often observed 
as jets, and it is generally believed that a toroidal mag- 
netic field is responsible for their coUimation (Benford 
1978, Blandford & Payne 1982). When a protosteUar wind 
strikes the ambient medium, a bipolar molecular outflow 
is produc ed. These outflows have been credited with the 
support (iNorman fc Silk 1980|. iMcKec 19891 iBcrtoldi & 



McKee 1996[) as well as the disruption (e.g.. Bally ct al 



1999D of molecular clouds and the dense clumps within 



them that produce star clusters. We seek a model of pro- 
tostellar outflows that is sufficiently detailed to permit a 
quantitative study of both their supportive and destruc- 
tive roles in the lives of their parent clouds. 

Bipolar molecular outflows display a number of common 
features that must be reproduced in any viable model. As 
summarized by Lada & Fich (1996), these include a nearly 
linear position- velocity relation (a "Hubble law") and a 
mass-velocity relation dM/dv cx with F ~ —1.8. It 
is debated whether these features result primarily from 
turbulent mixing in the region affected by the jet (e.g.. 



Raga et al. 1993| , [Stahlcr 1994| ) , or from the dy namics of 



a sho cked shell bounding the wind cocoon (e.g., Shu et al 
199li Masson & Chernin 1992, 1993). In this paper we 
investigate shells of ambient material set into motion by 
hydromagnetic protostellar winds, and demonstrate that 
these naturally produce the primary characteristics of ob- 
served outflows. (We do not address traces of high- velocity 
molecular gas sometimes found within these shells.) 

Our analysis of the wind force distribution (sQ follows 
Shu et al. (1995) and the suggestion made by Ostriker 
(1997); our investigation of shell motion (§||) follows Shu et 
al. (1991), Masson & Chernin (1993), and Li & Shu (1996). 
We generalize these results and show ([^), contrary to the 
conclusions of Masson & Chernin, that a combined model 
can reproduce the observed mass-velocity relationship as 
well as the position-velocity law. 



2. FORCE DISTRIBUTIONS FROM DISK WINDS AND 
X-WINDS 

We seek the distribution of wind momentum flux on 
scales of the dense clumps in molecular clouds, far larger 



than any scale associated with an accretion disk whose 
wind produces the outflow. The angular distribution of 
the wind momentum injection rate p^, at time t can be 
written as 
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where fi = cos 9 labels directions from the outflow axis. 
We wish to determine P{n), the normalized force distri- 
bution, which we assume is constant in time. 

Shu et al (1995) have shown that oc l/{rsm9)^ is a 
good approximation for X-winds. Since the wind velocity 
Wu, is approximately the same on different streamlines in 
this model, it follows that Pwv'^ oc l/(r sin O)"^. In fact, this 
should be a reasonably good approximation for any radial 
hydromagnetic wind that has expanded to a large distance, 
as the following heuristic argument indicates. Such winds 
expand more rapidly than the fast magnetosonic velocity 
Cf = B /(4:TTpi^y^^ , where we have assumed that thermal 
pressure is negligible. At large distances, the field wraps 
into a spiral with B ~ B^. First, consider the flow along 
streamlines. Flux conservation gives 2TrrArB^ = const in 
radial flow. Since is about constant at large distances, 
it follows that Ar ~ const, so that oc 1/r. Since 
Pw oc 1/r^ in a constant velocity, radial wind, it follows 
that v^/cf is approximately constant along a streamline 
at large distances. The flow approximates an isothermal 
wind, for which (w^/c/)^ increases logarithmically with 
distance from unity near the source. This implies that 
at large distances (uu,/c/)^ is approximately a constant 
across streamlines as well, so that Pwv'^ B^. But note 
that the field must become approximately force free at 
large distances: balancing the tension — _B|/(47rtn), where 
■CD = rsinO is the cylindrical radius, against the pressure 
gradient —{l/8Tr)dB^/dw gives B^ oc 1/w = l/rsm9. 

We conclude that PwVw oc i?^ oc l/(rsin0)^ is a general 
characteristic of radial hydromagnetic winds. 

This argument can be made more precise using a result 
due to Ostriker (1997). Suppose the wind arises from a 
Keplerian disk, where the wind density varies with initial 
radius vjq as po oc vUq'^ and the Alfven velocity varies with 
orbital velocity at the disk. The value q = 3/2 corresponds 
to the solution of Blandford & Payne (1982), whereas val- 
ues 0.5 < g < 1 were considered by Ostriker (1997). As- 
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sume that each streamUne expands to tu ^ tuo, so that 
the wind is significantly super- Alfvenic. The conserva- 
tion of specific energy, angular frequency, and mass flux 
along streamlines, along with the "isorotation" relation 

between B and v, give PwVw ^ C'{r)zulj'~'^^^'^w~^ , where 
C(r) is a slowly varying function of r (Ostriker 1997). We 
see that the heuristic argument above is valid provided 
that the disk radius njg, which varies from tuo.in to njQ^out, 
say, has a smaller range of variation than ro, which varies 
from an innermost radius r^corc to r. More precisely, if 
there is a power law relation between wq and sm9, then 



P{fj.) oc r^pu 
P{p) oc (sing) 



gives 



-2(l+£). 



SO we recover our earlier result if ^coro 
which allows the flow to be quasi-radial. 
ra7o,in = 1^0, out and the heuristic result 
accurate. Conditions at the axis set ^core- 
inner boundary is generally not considered in disk wind 
models, the Shu et al. (1995) theory posits a core of open 
field lines from the pole of the accreting star. The bal- 
ance of magnetic pressure and tension in the fiducial X- 
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ln(6'coro) 

(2) 

.out ; 

For an X-wind, 
should be quite 
Although the 



wind model gives Wc 



1.2 X 10' 



pc 



Hi 



2.5(1 + O.lSlogiorpc) AU, or 
0.18 log^Q Tpc), at a distance of 



parsecs. 



We have evaluated the accuracy of equation (Q) by cal- 
culating the wind force distribution analytically using the 
method outlined by Shu et al. (1995), as Ostriker (1997) 
suggested. We present this calculation in Matzner (1999), 
where we find that equation (^ is a good approximation 
at all angles more than fifteen degrees from the equator, 
and best matches the actual solution between one and ten 
degrees from the outfiow axis. The fiducial X-wind model 
(Shu et al. 1995) has P{ii) sin^ 9 greater by 40% toward 
the equator, for reasons not considered in equation (||); 
this corresponds to e ~ —1/50 within 10° of the axis at a 
distance of 0.1 pc. 

Although the approximation P{fi) oc (sin0)~^(^+'' is 
valid for ideal axisymmetric winds, we expect the force 
distribution to fiatten within some angle > ^coro in re- 
ality. This flattening could be produced by any of the 
mechanisms thought to create Herbig-Haro objects, e.g., 
jet precession, internal shocks from a fluctuating wind ve- 
locity, or the magnetic kink instability. Assuming that e is 
negligible and that ^ Ij we can therefore approximate 
the force distribution of a magnetized protostellar wind as 



\ii{2/9o){l + 9l-p^) 



(3) 



where the prefactor assures P{fi)dp — 1. Outflows may 
require a larger value of 6*0 than appropriate for winds 
themselves, if mixing between sectors (neglected here) di- 
lutes the momentum on axis; the current theory will still 
apply, with this larger Oq. However, so long as 0o ^ Ij 
the wind force is tightly concentrated along the axis: The 
formation of jets is an inevitable consequence of a hydro- 
magnetic wind. 



3. SHELLS DRIVEN BY PROTOSTELLAR WINDS 

We shall now explore the structure and motion of a shell 
of ambient material struck by a wind with the force distri- 
bution given by equation (^. Following Shu et al. (1991), 
Masson & Chernin (1992), and Li & Shu (1996), we ideal- 
ize the swept-up shell as thin and momentum-conserving. 
This is justified because both shocks are radiative for pro- 
tostellar wind velocities (Koo & McKee 1992a; however, 
see ^ for a consideration of magnetic pressure). We will 
also adopt the assumption that the fiow is entirely radial, 
so that mass and momentum are conserved in each angu- 
lar sector and there is no relative motion of the shocked 
fluids. 

A shell is driven by a "heavy" wind if less ambient ma- 
terial than wind material has been swept up; such shells 
travel at nearly the wind velocity, and the crossing time 
of the wind is therefore comparable to the outfiow age. 
Alternatively, a shell driven by a "light" wind is one that 
has swept up more ambient gas, and has decelerated signif- 
icantly. In the limit of a very light wind, both the wind's 
mass and its flight time can be neglected; this limit is ap- 
proached rapidly once a comparable mass has been swept 
up (Koo & McKee 1992b). Molecular outflows expand 
five to twenty times more slowly than their driving winds, 
and are comparably more massive. We may therefore ne- 
glect the wind's mass and flight time, and integrate the 
equation of momentum conservation in each direction 
dpw/dil = VsdMa{Rs, p) / dVt, where Ma is the ambient 
mass inside the shell radius Rg- We assume that the am- 
bient gas has a density pa — poir~^''Q{p), where poi is 
a constant and the angular factor Q{p) is normalized so 

that Jp Q{p,)dp — 1. We flnd that the shell radius is 



R 



^ (4-fcp)(3-fcp)P(M) 

^TTpOlQip) 



Pu,{t')dt'. (4) 



If the wind momentum is a power law in time, oc f'", 
then the shell velocity is given by v{p) — rjRg/t, where 

77 = 91ni?,(/i,t)/aint = (77in + l)/(4-A:p). (5) 

Equation (^) shows that the shell expands self-similarly, 
regardless of the wind history: its radial and velocity struc- 
tures are fixed, while its scale expands as J p^dt increases. 
Self-similarity is expected, because we have chosen a scale- 
free medium; other radial scales, such as the scale of the 
light-heavy wind transition, the wind coUimation scale, 
and the cooling length, are all small compared to a typical 
outfiow. 

Shu et al. (1991), Masson & Chernin (1992) and Li & 
Shu (1996) have all considered a steady wind {rjin = 1) 
and an ambient distribution appropriate for pre-stellar 
cores at the point of collapse: kp = 2, and Q{p) larger 
toward the equator, because of magnetic or rotational 
flattening. However, this is appropriate only in the re- 
gion whose gravity is dominated by the pre-stellar core, 
or where it prescribes a density higher than the ambient 



density: r ^ 0.07(crth/0.2 km s-^){nH/10'^ 
using the theory of Shu (1977). ^ 



pc. 



^Note that if the wind, star and core masses are each within about a factor of ten of the last, then the outflow must have emerged from its 
core if the wind is to be hght near its axis. This follows from the overwhelming factor (~ 10'') by which the axial force is enhanced; it is only 
exacerbated by any flattening of the core. The core mass may continue to affect the low-velocity, equatorial flow, however. 
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Outside this radius, anisotropics in the ambient gas are 
unHkely to correlate with the outflow direction, so we may 
assume ~ 1. We also then expect fcp ~ 0, 1 or ^ 2, 

if the lobe in question is smaller than, comparable to, or 
emerging from its parent star- forming "clump" . Assuming 
the ambient medium is isotropic and the outflow has not 
escaped its clump, and that e ^ 1 so that its efi^ect can be 
neglected in the outflow shape, our model reduces to: 



R 



'head 



= [i + (i-mW] 



2^fl-21-l/(4-fep) 




(6) 



where the radius of the lobe head, i?hcadj expands accord- 
ing to equations and (^) with ^ = 1. 

4. COMPARISON WITH OBSERVATIONS 

The shell described above is a Hubble flow in the sense 
that Vs{n,t) = riR,s{fJ.,t)/t. Since the relative hne-of- 
sight velocity is related to the line-of-sight distance by 
'^obs — fl^ios/t, the position-velocity (PV) diagram along 
the extent of an optically thin outflow is the same as its 
image to an observer situated in the plane of the sky along 
the short axis of the outflow. Whereas inclination causes 
a foreshortening of the outflow in a sky map, it causes the 
PV diagram to rotate. An elongated outflow (6*0 ■€i 1) will 
display a nearly linear PV diagram from the greatest to 
the least values of velocity and position. The agreement 
of self-similar outflows with the observed Hubble law was 
pointed out by Shu et al (1991) for the particular case 
kp — 2 and rji^ — 1; here we see that it is a quite general 
property. 

The mass- velocity relationship, dM/dvohs is a projec- 
tion of the PV diagram onto the velocity axis. Typical 
outflows show dM {vobs) / dvohs oc v^^^^, with F ~ —1.8. 
In an elongated outflow of inclination i, Vohs — vcosi, 
where v = |v|, because all but the lowest velocities are 
achieved at small angles from the outflow axis. Therefore, 
r = d\n{dM / dv) / d\n(v) except at low velocities. Because 
dM/dv = {dM/dij)/{dv/dn) oc w3-'=p/(^t,(5+<^-fcp)/(i+^)) 
from equations (|^) and (Q) with Q{fi) — 1, we find 



r 



' 1 + e ' 



for Dobs > 2w,„i„ where Umin ^ Oq 

imum space velocity (see eq. 0). Essentially all of the 
(non-equatorial) flow shows a value of F very close to —2 
(if e 1), in excellent agreement with a typical observed 
value of —1.8. The slope is shallower for Wobs <^ v^in'- be- 
cause dM/dvohs is a symmetric function of Vobs, F 
when Vohs = 0. The velocity at which dM/dv oc v~'^ 
fails is greater for more inclined outflows, or for larger 
^o; this could in principle constrain the inclination. The 
ability of our model to fit observational mass-velocity 
curves is demonstrated in Figure |l| for L1551, NGC2071, 
and NGC2264G. The qualitative agreement of X-winds in 
isothermal toroids with observed mass-velocity curves has 
previously been shown by Li & Shu (1996). Again, we 
see that this is a general feature of momentum-conserving 
shells driven by hydromagnetic winds, and is not linked 
to particular models for the source of the wind nor the 
ambient medium, so long as this medium has a power law 
density distribution. 



2/(4-fcp) 



(7) 



VhQud is the min- 



The broadening angle can be constrained using the 
range of velocities for which F ~ — 2: this law holds for 
at least a factor of 4 in the outflows shown in Figure 1. 

2/f4 k ) 

This must be less than about Whead/(2wmin) = '' ' 

so 6q <^ l0^i-^(i~'^p/4). However, the velocity factor is also 
reduced by inclination at a given Oq; in Figure 1 we show 
that 6*0 = 10~^ allows an inclination of ^ 45° to fit the 
data; 6q = 10^^ '^ gives too small a velocity range unless 
i = 0. 

If Pobs, lobe is the net momentum of an outflow lobe along 
the line of sight, another constraint on 9o comes from the 
ratio. 



Pobs, lobe 



^^^dM/dvohs) head 



21n(go-^) 
4-fc„ ' 



(8) 



valid if e, ^0 ^ 1- The result was obtained by integrating 
the momentum along the line of sight for an uninclined 
(i = 0) outflow; we did not assume that /i ~ 1, but in- 
stead used the exact expression Wobs = /if appropriate 
for each shell. It is interesting to note that this agrees 
exactly with the expression ln(uhGad/wmin) one would es- 
timate from dm/dv oc with ^ — \. Although the 
ratio was derived for zero inclination, it is actually in- 
dependent of I, because the numerator and denominator 
scale together. Observations place a lower limit on the 
ratio (an upper limit on ^o), since the net momentum 
may be underestimated. We find that 0o < 10~^'^ for 



L1551 (iMoriarty-Schieven fc SncU 19881) and gp < lO"^- ^ 
for NGC2071 ( |Moriarty-Schievcn, Hughes fc SneU 1989| ), 
assuming kp — 0. Again, note that Oq 2± 10~^ is consistent 
with the data. 

The current model reproduces the typical extents an_d 
velocities of observed outflows. For instance, equations 
and (^) imply that a wind of momentum 20 Mq km s" 
with 9q — 10^^, blowing steadily into a uniform density 
10^ cm~3 for 10^ years, drives lobes whose heads decel- 
erate to 7.4 km s""'^ and expand to 1.5 pc (each) in this 
time. 

5. CONCLUSIONS 

Perhaps the most remarkable properties of bipolar 
molecular outflows, apart from their intensity and fre- 
quency in regions of active star formation, are their high 
degrees of collimation and the commonality of the rela- 
tion dM/dv oc v~^-^. We have shown that hydromag- 
netic winds are naturally coUimated so that the force 
distribution, PwV^^ cjc l/sin^0, approximately leads to 



dM/dv 



in any power-law ambient medium, pro- 



vided the interaction is momentum conserving. Our re- 
sults show that that the conclusions reached by Shu et al 
(1991) and (1995) and Li fc Shu (1996) for steady X-winds 
in media with 1/r^ density distributions are far more gen- 
eral, and in addition are in good agreement with observa- 
tions of protostellar outflows. 

The fact that outflows are often observed to have dM/dv 
slightly shallower than indicates that the theory is 
only approximate. A number of effects we have not con- 
sidered could lead to a deviation from the -2 slope: CO 
self-absorption at lower velocities, mixing of radial mo- 
mentum between angles, or the generation of lateral mo- 
mentum when the wind impacts the shell (Masson and 



10 
( km ; 



Fig. 



Mass- velocity relations for Ll^ 



qriarty-Schieven fc Snell 1988 ), NGC2071 (squares: Moriarty-Schieven, Hughes & 
Snell 1|89|), and NGC2264G (diamonds: |Lada fc Fich 1996| ), fit with an outflow model (eg. M: kp = 1, 9o = 0.01, e = 0, solid lines). The 
inclination has been chosen to make the spectral breaJj near the lowest observed velocities: 50° tor L1551, 40° for NGC2071 and NGC2264G. 
Filled symbols: blue lobes; open symbols: red lobes. Velocities are relative to vq, a center velocity, chosen to maximize the symmetry between 
each pair of lobes. Curves have been offset vertically for clarity. We have taken e = so that the slope F is —2 in the model; but see 3a. 



Chernin 1993). In the current model, F can differ from —2 
either because of the wind force does not exactly follow 
PwVw (sin0)~^ (e 7^ 0), or because the ambient medium 
is not a single power law {kp varies). The model also pre- 
dicts a flattening of dM/dv at low velocities, which could 
raise the estimate of F. 

Let us consider the possibility that e is to blame for 
F 7^ —2. From equation (^, F ~ —1.8 requires ~ 
5(4 — kp): the wind is slightly more concentrated toward 
the axis. Equation implies that the number of decades 
of disk radius required to give this value of e is approxi- 
mately logio(l/6'corc)/[5(g - 1)(1 - fcp/4)]. The Blandford 
& Payne (1982) model has g = 3/2; it would therefore re- 
quire about 2-4 decades of disk radius to give F = —1.8 



for 0corc ^ 10-5 



and 2 > fcp > 0. A disk with g < 1 



(initial wind density increasing outward, e.g., Ostriker's 
[1997] models), has its wind force weighted toward the 



equator relative to Pwvfu oc (sin( 



and produces out- 



flows steeper than dM/dv (x v~'^ . X-wind models share 
this behavior, as they predict e ~ —1/50. 

Our model assumes that the shocked wind and ambient 
gas form a thin shell. Although this i s appropriate for un- 
magnetized gas ( Koo fc McKee 1992a| ), the fact that winds 
are collimated magnetically raises the possibility that out- 



flows might become inflated with a cocoon of magnetically- 
supported shocked wind, before or after the wind shuts off. 
This depends on the wind's terminal Alfven Mach num- 
ber (inversely proportional to Poynting flux), and also on 
whether the field remains ordered or becomes tangled. Be- 
cause the Poynting flux decreases as the wind collimates 
(^corc decreases), and also as the kink instability develops 
(Choudhuri & Konigl 1986), it is reasonable to ignore the 
magnetic pressure of the shocked wind. 
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